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Abstract. - We study a random walk model in which the jumping probability to a site is 
dependent on the number of previous visits to the site, as a model of the mobility with memory. To 
this end we introduce two parameters called the memory parameter a and the impulse parameter p. 
From extensive numerical simulations, we found that various limited mobility patterns such as sub- 
diffusion, trapping, and logarithmic diffusion could be observed. By the memory, a long-ranged 
directional anti-correlation kinetically-induces anomalous sub-diffusive and trapping behaviors, 
and transition between them. With random jumps by the impulse parameter, a trapped walker 
can escape from the trap very slowly, resulting in an ultraslow logarithmic diffusive behavior. Our 
results suggest that the memory of walker's has-beens can be one mechanism explaining many of 
empirical characteristics of the mobility of animated objects. 



q Introduction. — Since Einstein's study on the mo- 
•^tions of small particles suspended in liquids [l], diffusion 
£>Land random walk have been a paradigm in describing and 
modeling the mobility involving seemingly irregular mo- 
t CU ions exhibited by both material and animated objects [2]. 

A multitude of variations of the random walk model with 
t-H additional details have been widely studied [3j[4], Gover- 
ning diverse physical and biological domains ranging from 
charge transport in disordered conductors to foraging pat- 



terns of animals and fishes. 
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^T*) From the random- walk perspective, even the movement 
of a human individual may look irregular and random. In 
contrast to material objects, however, humans are con- 
{S| scious beings and driven by motivations by their will and 
purposes modulated by external cause, so their mobility 
patterns cannot be completely random. Indeed, recent ad- 
jvances in digital technology provided us with a plenty of 
i empirical data for human mobility 5,6, which revealed a 
number of key empirical characteristics of human mobil- 
ity patterns, three of which we highlight here: First, even 
though the individual trajectory bears some degree of ran- 
domness, most part of it is highly predictable, being em- 
bedded in a well-defined region in space with the radius of 
gyration of the trajectory growing logarithmically in time 
pi. Second, the mobile object (a human individual in this 
case) spends disproportionate times in different locations 
[7]. Third, at the population level, there is high degree of 
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heterogeneity across individual's degree of mobility [6]. 

In this work, our key premise is that at the heart of these 
human mobility patterns does underlie the role of memory: 
Conscious beings do have memory of where he or she has 
been and tend to revisit some of the have-beens, his or her 
favorite spots. By incorporating such a memory effect in 
the dynamic rule of a random walk model, we investigate 
the consequential role of memory in the mobility patterns. 

We will demonstrate that the memory induces a long- 
ranged directional anti-correlation, which kinetically slows 
down the walker's mobility, resulting in sub-diffusion and 
trapping, and transition between them. Logarithmic dif- 
fusion can occur by introducing the impulsive random 
jumps. Therefore, many of characteristics of limited mo- 
bility can be explained by simple rules based on memory. 

Model. — Our model is defined on a hypercubic lat- 
tice in d-dimensions. Starting from the origin, the walker 
takes a random walk with jumping probability depending 
on the number of previous visitations to the target sites. 
Specifically, we denote rii(t) the number of visitations to 
the site i upto t steps by a random walker starting from 
the origin at time t = 0. The jumping probability Wj^(t) 
from the site i to site j at time t is given by 
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with the normalization Z = J2/[ n j(t) + M a by the re ~ 
stricted sum over nearest neighbors j of the site z, and 
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Fig. 1: Typical trajectories from the model, corresponding to 
the parameter combinations (a,p) = (0.4,0) (a), (0.8,0) (b), 
and (1,0.2) (c), each representing the sub- diffusive, trapped, 
logarithmic diffusive behaviors, respectively. Trajectories upto 
t = 2 15 in two dimensions are shown, colorcoded by the visi- 
tation numbers (see the color-bar). Note the difference in the 
ranges of view field, as well as the range of variations in visi- 
tation numbers in the three cases. 



d the spatial dimension. Thus the jumping probability 
changes dynamically, depending on the history of the ran- 
dom walk. 

The two parameters a and p characterize our model, a 
determines the degree of memory-dependent revisitation 
tendency. So we call it the memory parameter. The higher 
a is, the walker tends to go back to sites previously visited 
more often. The parameter p is introduced to implement 
occasional purely random jumps, which may occur due to, 
e.g., impulsive decisions of a mobile walker. So we call p 
the impulse parameter. With p = 0, the walker always 
performs the memory-dependent move. With p = 1, the 
model is reduced to simple symmetric random walk. In 
Fig. 1, typical trajectories from our model in two dimen- 
sions with different combinations of (a,p)-parameters are 
shown. 

Random walk models in similar vein to ours had been 
studied in relation to self-interacting random walks 8- 
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For example, Duxbury and de Queiroz 10 assigned 
to each (simple) random walk trajectory the "energy" in 
the form of exp (— g J2i n i a ): where ni is the visitation 



numbers of the site i as in our model. When g > and 
a > 0, constrained trajectories are favored (note that even 
a = case is not reduced to simple random walk in this 
model). Depending on g and a, different behaviors occur, 
such as anomalous sub-diffusion, trapping at origin, or 
self- avoiding walks 
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Although apparently similar, this 
class of models are equilibrium- type models, whereas our 
model is a generative model with explicit dynamic rules. 

More recently, human mobility models directly moti- 
vated by data from digital records of human travels and 
mobility have been proposed |13]jl5] . For example, Song 
et al. derived two rules of human mobility from the mobile 
phone usage data 15 . They introduced a human mobil- 
ity model incorporating those rules, which they termed 
the preferential return and the exploration, similar to our 
memory and impulse parameters. With such data-driven 
modeling, they explained many characterisics of empirical 
human mobility patterns, e.g., the ultraslow logarithmic 
growth of the radius of gyration of individual's trajectory 
in time. This model and ours share key properties in dy- 
namic rules, but the precise implementation is somewhat 
different. First, the model of Ref. [15] is a spatial network 
model, whereas ours is a lattice model. Second, the pa- 
rameters in our model cover a broader range, therefore we 
can address the generic role of these parameters, which 
might be relevant for novel mobility patterns associated 
with memory effect, other than those identified thus far. 

Role of memory parameter. — The memory pa- 
rameter a determines the degree of memory-dependent 
revisitation tendency. We are primarily interested in the 
cases of positive a, for which the walker tends to revisit 
the sites previously visited many times, resulting in a con- 
strained trajectory [see Fig. 1(a) for example] compared to 
simple symmetric random walks, similarly to interacting 
random walks [9llTo] . When a = 0, the model becomes 
memoryless and is reduced to simple symmetric random 
walk. When a is negative, the walker disfavors revisita- 
tion, which would lead to a stretched trajectory, similarly 
to weakly self- avoiding walks M. In the limit a —> — oo, 
the model becomes the self-avoiding walk. In the oppo- 
site limit a — >• oo, which is of more interest to us, the 
walk would become so constrained that the walker will be 
trapped at the origin [see Fig. 1(b) for example]. The in- 
teresting questions are the existence of phase transition to 
the trapped trajectory at finite a c and how the normal 
diffusive behavior at a = and trapping at a = a c are 
interpolated as a increases. Thus in the following we focus 
on positive a regime. 

To answer these questions, we perform extensive nu- 
merical simulations. To focus on the role of memory pa- 
rameter, we fix p = and vary a in the range from 
to 1 (Fig. 2). We measured the average dispersion of the 
walker in terms of the mean square displacement, (r 2 (t)), 
as a function of the number of steps t. (• • • ) denotes the 
ensemble average. The dispersion increases with number 
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Fig. 2: Numerical simulation results for < a < 1 and p = in 
one dimension, (a) Mean square displacement (r 2 (t)) as a func- 
tion of number of steps t are displayed in double-logarithmic 
scale. Straight lines are guidelines with slopes 1 (top) and 0.5 
(bottom). Different symbols correspond to different a from 
(top-most) to 1 (bottom-most) in the interval of 0.1, from top 
to bottom. (Inset) The exponent v vs. the memory parame- 
ter a, obtained by fitting the data to Eq. (2) for large t. (b) 
Double-logarithmic plot of the probability of visiting an unvis- 
ited site at step £, (A(£)), as a function of t. Straight lines are 
guidelines with slopes —0.5, —1, and —1.3, from top to bottom. 
Same symbols as in (a) are used. (Inset) The exponent a vs. 
a, obtained by fitting the data to Eq. (3) for large t. Ensemble 
averages were taken over 10 4 independent trajectories. 



of steps in a power-law manner at long times, 



{r\t))~t 



2v 



(2) 



The exponent v is 1/2 for ordinary diffusion, and v < 1/2 
indicates sub-diffusive behavior. We found that in one di- 
mension, v decreases continuously from 1/2 at a = to 
zero near a c ~ 0.5 [Fig. 2(a)]. This suggests a trapping 
transition at finite a c ~ 0.5 and sub-diffusion with contin- 
uously varying exponent v for < a < a c . 

We also measured the average probability that the 
walker visits to an unvisited site at step t, (A(t)). For 
one dimensional simple symmetric random walk (a = 0), 
it is known [3] that (A(t)) decreases with t as a power law, 



<A(t)> ~ t~ 



(3) 



with a = 1/2. We found that a increases with a, be- 
comes a = 1 at a = a r ~ 0.5, and further increases for 



a > a c [Fig. 2(b)]. This indicates that the average total 
number of visited sites (S(t)) upto t- steps, given by the 
sum of (A(t)) over t, exhibits distinct behaviors depend- 
ing on a. (S(t)) diverges as ~ ^-^O) f or a < ac (diffus- 
ing), whereas it becomes finite for a > a c (trapped). At 
the critical point a = a c , (S(t)) increases logarithmically. 
This result is also in support of the existence of a trapping 
transition at finite a c . 

Thus, with the memory parameter we can obtain the 
slow or limited mobility in walker's motion. Depending 
on its strength, it can kinetically induce anomalous sub- 
diffusion or trapping, and transition between them. The 
heterogeneity of individual trajectory would follow at the 
critical point. It can also be thought that the memory 
parameter a can take different values person to person, as 
it characterizes an individual's degree of revisitation ten- 
dency. However, with the memory parameter alone, the 
walker is either diffusing or gets trapped at some localized 
region in space. Ultraslow dispersion can only be expected 
for a single specfic value of the memory parameter at the 
critical point. 

Role of impulse parameter. — The impulse param- 
eter p models the degree of impulsive, random decisions 
in walker's motion, introducing occasional unbiased ran- 
dom jumps. Intuitively, introduction of unbiased random 
jumps (p > 0) will give rise to larger dispersion of the 
walker's trajectory than p — case. A question of interest 
is what the effect of random jumps is to a trapped walker. 
To answer this question, we fix the memory parameter to 
be a = 1, high enough for a strong trapping, as well as 
in accordance with the empirical data 15 , and vary the 



impulse parameter p in the range from to 1 (Fig. 3). 

By measuring the mean square displacement (r 2 (t)) as 
a function of t, we find that there exists a critical impulse 
parameter p c w 0.5, above which the walker is no longer 
confined in the trap and diffuses [Fig. 3(a)]. For p > p c , 
v increases continuously from zero and becomes 1/2 for 
p — >• 1, crossing-over to normal diffusion. This result is 
corroborated by the average probability of visiting an un- 
visited site (A(t)) decaying as Eq. (3) with the exponent 
<j varying from 1 to 1/2 as p increases from p c w 0.5 to 1 
[Fig. 3(b)]. 

The role of impulse parameter for p < p c is more intrigu- 
ing. Meausuring ^, we have v ~ for p < p Cl suggesting 
that the walker is still trapped [Fig. 3(a)]. However, a 
is stuck at a = 1, that is, (A(t)) ~ t -1 , regardless of 
p < p c [Fig. 3(b)]. This means that the average total 
number of visited sites (S(t)) still increases logarithmi- 
cally, (S(t)) ~ logt, suggesting that the walker is in fact 
marginally escaping from the trap. Such behavior can be 
seen from the example trajectory shown in Fig. 1(c), where 
one can identify more than one traps. 

To illustrate mobility patterns exhibited by different 
phases, we calculated the probability distribution P;(£), 
defined as the average probability to find the walker start- 
ing from origin at site i after t-steps (Fig. 4). For the 
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Fig. 3: Numerical simulation results for a = 1 and < p < 
1 in one dimension, (a) (r 2 (t)) vs. t in double-logarithmic 
scale. Straight lines are guidelines with slopes 1 (top) and 0.5 
(bottom) . Different symbols correspond to different value of p 
from 1 (top-most) to (bottom-most) with the interval of 0.1, 
from top to bottom. (Inset) The exponent v vs. a, obtained by 
fitting the data to Eq. (2) for large t. (b) Double-logarithmic 
plot of (A(t)) vs. t. Straight lines are guidelines with slopes 
0.5, 1, and 1.3, from top to bottom. Same symbols as in (a) 
are used. (Inset) The exponent a vs. a, from fitting the data 
to Eq. (3) for large t. 

memory-only model (p = 0), we have diffusing (a < a c ), 
critical (a = a c ), and trapped phases (a > a c ) [Figs. 4(a- 
c), respectively]. In particular, in the trapped phase, Pi(t) 
becomes time-independent at long times, suggesting the 
walker is competely trapped, visiting a set of trapping sites 
repeatedly ad infinitum. With the impulse parameter, the 
behavior is similar for the diffusing phase [Fig. 4(d)]. For 
P <Pa the random jumps by the impulse parameter leads 
to larger dispersion in Pi(t) [Figs. 4(e,f)]. In particular, 
Pi(t) is not strictly stationary for < p < p c but gets 
dispersed with time very slowly (in logarithmic orders) 
[Fig. 4(f)]. 

Such logarithmic dispersion is reminiscent of the Sinai 
diffusion in quenched-disordered media [18]. The revis- 
itation tendency kinetically generates random traps, at 
which the walker is trapped in short timescales; thanks 
to random jumps, however, the walker can slowly escape 
from the trap, the timescale of which is of exponential or- 
der in the trap strength [l8] , but the walker will soon get 
re-trapped at another kinetically-generated trap. Such a 
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Fig. 4: Color-coded plot of the time evolution of the probabil- 
ity distribution Pi(t) that the walker will be at site i at step 
t, for parameter combinations (a,p) as (a) (0.2,0), (b) (0.5,0), 
(c) (0.8,0), (d) (0.8, 1), (e) (0.5, 1), and (f) (0.2, 1), representa- 
tive of the sub-diffusive (a, d), critical (b, e), trapped (c), and 
logarithmc diffusive (f) regimes. Note the logarithmic scale in 
time directions. 

series of trapping leads to the observed ultraslow disper- 
sion. In our model the trap strength changes dynamically, 
making it even slower than ~ (logt) 4 dispersion for the 
original Sinai diffusion with static disorder [18] . 

Kinetic origins of the anomalous diffusion. — 

Our numerical simulation results show that the memory 
parameter a can modulate the diffusivity of the mobility 
with memory, from sub-diffusion to trapping. In order 
to gain more microscopic understanding how such limited 
mobility is generated by the memory, we measured the di- 
rectional correlation function, C(s), denned as the average 
correlation between the directions of two steps separated 
by the lag of s steps, that is, C(s) = (rj(t)r)(t + s)), where 
rj(t) = r(t) — r{t— 1) is the step direction in one dimension 
and the overbar denotes the time average. In simple sym- 
metric random walk (a = 0), C(s) = 0. We found that 
for a > 0, this directional correlation is non-zero. 

In the sub-diffusive regime (0 < a < 0.5), two notable 
features emerge for C(s) [Fig. 5(a)]. First, C(s) is on 
average negative in the entire range of s, and its abso- 
lute value for small s increases with a. This means that 
the walk becomes more and more strongly anti-biased to- 
wards the current position anytime. Second, |C(s)| de- 
cays as a power law, \C(s)\ ~ with a crossover be- 
tween two regimes. The power-law exponent is ( ~ 1/2 
for small s and £ w 2 for large 5, and the crossover 
timescale decreases with a. This long-ranged directional 
anti-correlations might be the microscopic origin for the 
anomalous sub-diffusive behavior, with a similar mecha- 
nism as in the continuous-time random walks with anti- 
correlated fractional Gaussian noise [17]. Meanwhile, in 
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Fig. 5: Detailed kinetics of the model, (a) Directional corre- 
lation function C(s) as a function of the step lag s in double- 
logarithmic plot. Straight lines are guidelines with slopes —0.5 
(dotted) and —2 (dashed), (b) Block waiting time distribution 
?Pb(t) with the blocksize £ = 7 in double-logarithmic scale. 
Straight line is a guideline with slopes —1. Same set of sym- 
bols are used in (a) and (b). 

the trapped regime (a > 0.5), the directional correlation 
develops an alternating behavior around zero, reflecting 
the back- and- forth motion within the trap [Fig. 5(a), in- 
set]. The amplitude of alternation increases with a, as the 
trapping becomes stronger. 

A well-known mechanism for sub-diffusion is the power- 
law distribution of waiting times between consecutive 
steps, or pausing times [19]. When the unbiased ran- 
dom walk with finite step length has the waiting time 
distribution ip(r) of a power-law form, tp(r) ~ r ~^ 1+/i \ 
with < \i < 1, its dispersion becomes sub-diffusive as 
(r 2 ) ~ t M at long times. To consider how such an effect 
manifests in our model, we divide the lattice into non- 
overlapping linear blocks of size t (typically t = 7, to cover 
the typical size of a trap), and calculated the histogram 
of the time interval r between two consecutive crossings 
of the block boundaries by the walker, which we call the 
block waiting time distribution i/jb(t). 

We found that in the trapped regime (a > 0.5), 
the block waiting time distribution follows a power law, 
^b(t) ~ t -1 , for large r [Fig. 5(b)], leading diverging av- 
erage waiting time. This divergence is likely rooted from 
the existence of strong trapping, with diverging mean es- 
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Fig. 6: Color-coded plots of the exponent v in (a,p)-parameter 
space in (a) one and (b) two dimensions. Dashed lines are 
guidelines of a(l — p) = 1/2. 

caping time. The power-law behavior with [i — in the 
trapped regime is apparently consistent with the lack of 
dispersion v = in that regime. In the sub-diffusive 
regime (0 < a < 0.5), however, the block waiting time 
distribution is found to decay exponentially, rather than 
as a power law [Fig. 5(b)]. This result shows that the sub- 
diffusive mobility of the walker in our model is quite differ- 
ent from that of an unbiased random walker with power- 
law waiting time distribution; our walker is not waiting at 
particular locations but moves more actively from block 
to block constantly. It is the anti-correlated motion re- 
vealed above that provides key kinetic mechanism behind 
the limited mobility. 

In two-dimensions. — We performed numerical sim- 
ulations of our model in one dimension for other combi- 
nations of (a, p)-parameters. We found that the critical 
value of memory parameter a c depends on the impulse 
parameter p; starting from a c « 0.5 at p = 0, a c increases 
with p [Fig. 6(a)]. We also performed numerical simu- 
lations in two dimensional square lattice, finding overall 
similar behaviors. A notable difference is that the log- 
arithmic dispersion for a > a c and p > is seen more 
prominent in two dimensions [Fig. 6(b)]. Our numerical 
data suggest the dispersion in the logarithmic regimes as 
~ (logt) 7 with 7 w 0.3 in one dimension and 7 ~ 1.5 in 
two dimensions. This result shows that the effect of ran- 
dom jumps by the impulse parameter becomes manifestly 
more significant in higher dimensions. 

Continuum limit. — The discrete-time master equa- 
tion for Pi(t) can be written as, 

Pi(t + 1) - Pi(t) = J2' Mt) p i(t) - v>i,i(t)Pi(t)] . ( 4 ) 

3 

where the transition probability Wj^(t) is given by Eq. (1) 
and the restricted summation runs over nearest neighbors 
j of i. Decomposing Wjj into memory-dependent and ran- 
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dom jumps, we can rewrite the above equation as 

Pi(t + 1) - Pi(t) = (1 - P )J2' ImAQPjV) ~ hMpM 



2d 



(5) 



Here fi jti (t) = MQ/EkM*) with M*) = M*) + 1]", 
denotes the memory-dependent jumping probability. Let 
us consider a continuous time and space limit of Eq. (5). 
One can rewrite ^j^(t) in one dimension as = 
- ± - etc., assuming 



l 

2 2 ~ 2 2 arc 

that t), the continuum analog of is analytic at 

This leads to the continuum limit of Eq. (5) in one 



dimension as 
dP(x,t) _ ld 2 P(x,t) 



dt 



dx 2 



(i-p) 



d_ 

dx 



P(x,t) 



d(f){x,t)/dx 



<j>(x,t) 

Taking advantage of the relation between <f>(x, t) and 
P(x,t), one finally obtains a highly-nonlinear integro- 
differential equation for P(x,i), 



dP(x,i) 



dt 



ld 2 P(x,t) 
2 



dx 2 



a(l-p) 



d_ 

dx 



P ^ X ^dx~ X ° g 



P(x,t')dt' | . 



(7) 



Eq. (7) is too nonlinear to be solved exactly, yet it can 
be used to locate the transition point between diffusion 
and trapping as follows. Approaching to the transition 
point, P(x, t) starts to become time-independent. Looking 
for the condition of time-independent P(x,i) for large t, 
equivalently n(x,t) « tP(x), we find the transition point 
to trapping as a c (l — p) = 1/2, which is in good agreement 
with the numerical simulation results (Fig. 6). 

Summary and discussion. — We have studied a ran- 
dom walk model in which the jumping probability to a site 
is dependent on the number of previous visits to the site, 
as a model of the mobility with memory, such as that 
of a human individual or an animal. Through extensive 
numerical simulations, we found that by varying two pa- 
rameters called the memory parameter a and the impulse 
parameter p, various limited mobility patterns such as the 
sub-diffusion, trapping, and logarithmic diffusion could be 
observed. With the memory parameter a, a long-ranged 
directional anti-correlation is developed in the walker's 
motion, resulting in anomalous sub-diffusive behaviors. At 
the critical value of a c , the transition from sub-diffusive to 
trapped motion occurs, above which the walker's trajec- 
tory becomes confined in a localized trapping region. By 
allowing random jumps with the impulse parameter p, the 
walker can escape from the trap, albeit very slowly, but 
gets trapped to another trap kinetically-generated by the 
walker's own motion, resulting in an ultraslow logarithmic 
diffusive behavior. 



Our results demonstrate that the memory of walker's 
has-beens can be a mechanism explaining many of empir- 
ical characteristic patterns of human mobility [5j|7] . The 
revisitation tendency driven by the memory has also been 
proposed and supported by empirical data in the more de- 
tailed, data-driven modeling of human mobility 15 . Our 



work extends this approach to encompass broader spec- 
trum of limited mobility patterns, including sub-diffusion 
and trapping, and the transition between them. It would 
be interesting to see if our modeling framework with the 
memory and impulse parameters could be applied to vari- 
ous empirical data for human and animal's mobility, which 
might uncover further universality and specificity in the 
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